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While detection of the "local form" bispectrum of primordial perturbations would rule out all 
single-field inflation models, multi-field models would still be allowed. We show that multi-field 
models described by the SN formalism obey an inequality between /nl and one of the local-form 
trispectrum amplitudes, tnl, such that tnl > f (§/nl) 2 with a possible logarithmic scale depen- 
dence, provided that 2-loop terms are small. Detection of a violation of this inequality would rule 
out most of multi-field models, challenging inflation as a mechanism for generating the primoridal 
perturbations. 

PACS numbers: Valid PACS appear here 



Can we rule out inflation as a mechanism for gener- 
ating primordial curvature perturbations? Inflation is 
indispensable for explaining homogeneity and flatness of 
the observable universe Yet, its predictions for the 
statistical properties of primordial curvature perturba- 
tions may be falsifiable. 

The basic predictions that inflation generates adia- 
batic, nearly scale-invariant, and nearly Gaussian pri- 
mordial curvature perturbations 0, Q are all consis- 
tent with the current observations (see, e.g., jH). No- 
tably, many inflation models predict that the amplitude 
of fluctuations on large scales is greater than that on 
small scales. In terms of the power spectrum of pri- 
mordial curvature perturbations £, we say k 3 P^(k) oc 
fc" s_1 with n s < 1. The power spectrum is defined by 
(CkCk') = (2n) 3 S(k + k')Ps(k). The latest observations 
give n s = 0.96 ± 0.01 0, Hf, which may be taken as evi- 
dence for inflation. 

The future, more sensitive experiments will continue to 
test the other predictions: adiabaticity and Gaussianity 
of fluctuations. In this paper, we shall focus on the lat- 
ter. Departure from Gaussianity, called non- Gaussianity, 
has emerged as a powerful test of inflation over the last 
decade (see [f| for reviews). 

One of the major theoretical discoveries made from 
these studies is that all single-field inflation models yield 
a specific amount of non-Gaussianity in the so-called 
squeezed limit of the bispectrum (Fourier transform of 
the three-point correlation function) of C, given by /nl = 
jTj(1 — n s ) ~ 0.02 (also see Q). Here, /nl character- 
izes the amplitude of the so-called "local form" bispec- 
trum UGH: 



B C = ^/nl [P c (fci)P c (*2) + (2 perm.)] , (1) 



where OlLiC^)) = (27r) 3 <5 3 (y\ k t )B c (k u fc 2 , k 3 ), and 
the "squeezed limit" is given by taking ^3 <C k\ » k 2 , 
i.e., B ( (k 1 ,k 2 ,k 3 ) -> fhLP((ki)Pdh)- All single-field 
inflation models predict (1— n s )P^(ki)PQ (k 3 ) in this limit. 

The current best limit is /nl = 32 ± 21 (68% CL; 0]). 
As various second-order effects generate /nl = 0(1) (see 
11] for a review and references therein), a convincing de- 
tection of /nl ^> 1 would rule out all single-field inflation 
models. The Planck satellite is expected to reduce the 
error bar by a factor of four (lo| . 

However, detection of /nl would not rule out multi- 
field models. How can we test them also? Our work in 
this paper is motivated by the Suyama-Yamaguchi in- 



equality, t N l > (|/nl) 2 [12j- Here, t N l is one of the 
amplitudes of the local-form trispectrum defined by [l3T ] 



T c =T NL [Pc(l k i 



)P f (k 3 )P c (ki) + (11 perm.)], (2) 



where (nti CO*)) =l^) 3 S 3 (E l Wc^i, k 2 , k 3 , fc 4 ). 

As emphasized in [llj, if the new experimental data 
(such as Planck) detect /nl (hence ruling out single-field 
models) but do not see tnl large enough to satisfy the 
above inequality, then a large class of multi-field models 
may be ruled out. The crucial question is then, "how 
generic is the Suyama-Yamaguchi inequality?" It was 
pointed out in ll| that this inequality may not be generic 
enough, as there are cases where this inequality is not 
satisfied. Recently, Suyama et al. 1 li considered the 
same issue, where they have truncated the SN expansion 
(given below) at the second order and have considered 
a part of 1-loop corrections. The goal of this paper is 
to find a more general inequality than theirs. We shall 
retain the terms up to the fourth order of SN expansion, 
as these terms are required for the consistent calculations 
up to the 1-loop level. As a result, we find a weaker bound 
than the original Suyama-Yamaguchi inequality. This is 



2 



relevant because, as shown in 15[, large and observable 
primordial non-Gaussianity can be generated when the 
loop contributions dominate over the tree contributions 
in the bispectrum and/or in the trispectrum. 

Throughout this paper, we shall consider a class of 
multi-field models which satisfy the following conditions: 

1. Scalar fields are responsible for generating curva- 
ture perturbations; thus, potential contributions 
from vector fields (see [Iff for a review and ref- 
erences therein) are ignored. 

2. Fluctuations in scalar fields at the horizon crossing 
are scale invariant and Gaussian. 

Therefore, we assume that non-Gaussianity is generated 
only on su per horizon scales, according to the SN for- 
malism [ii. Il7j. While the "quasi-single- field inflation" 
model proposed by Chen and Wang [18| is an example to 
which this condition may not apply, their model yields 
Tnl 3> /nl, satisfying the inequality. Yet, the condition 
2 is probably too strong. Whether this condition can be 
relaxed significantly merits further investigation. 

According to the SN formalism, the curvature pertur- 
bation, C, is |iven by derivatives of the number of e-fold, 
N(t,t*) = L Hdt\ with respect to scalar fields, ip a , at 
the horizon-crossing time (a*H* = k): 

C(x,i) = JV a (t,*0^(x) + iAU(M*)<^*(x)V*(x) • • • 

(3) 

where Sipl is a fluctuation of ip a evaluated at i*, i.e., 
<V£(x) = <V(i*,x). Note that N a = dN/dcp% and 
N ab = 8PN/d<f%d(£. 

The second condition above implies that the power 
spectrum of scalar fields is given by 



shall include all of the 1-loop contributions, while some 
of the higher-order loop contributions are also included. 
The power spectrum is given, up to the 4th order, by 

Vc, =V* \N a N a + Tr(N 2 )V* ln(fcL) + N a N abb V* ln(fc max L) 



+ ^N acc N abb V 2 ln(k mllx L) 
+ N abcc N ab P 2 ln(fcL) ln(fc max L) 



(5) 



where we have used the following notations: N a N a = 
J2 a N 2 a and Tr(7V 2 ) = £ Qb N ab N ab . The L is a finite size 
of a box which is chosen to be much larger than the region 
of interest, such that the condition Lk ^ 1 is satisfied 
for arbitrary k, and A: max is the ultra-violet cutoff. The 
1st term is the tree contribution; the 2nd and 3rd terms 
are the 1-loop contributions; and the 4th and 5th terms 
are the 2-loop contributions. 

This result can be simplified by using the following 
quantities (see Eq. (25) of [HI): 



N a = N a + ]-N abb V* ln(fc max L), 



(6) 



N ab = N ab + -N abcc V* ln(fc max L). (7) 

Then Eq. flSJ becomes 

V ( = N a N a P, (1 + 7> loop + ...), (8) 
where we have defined a positive-definite quantity 



Tr(A> 2 ) 

Pioop = ~ - L V*HkL). (9) 
N a N a 



(6<pl(U)8ipl,(Q) = (27r) 3 ^ 3 )(k + k')<5 ah ^^, (4) 

where V* = {H*/2tt) 2 . Note that we have assumed that 
scalar field fluctuations with different indices are uncor- 
rclatcd, {5ip a 6cp b ) oc S ab . This can be done without loss 
of generality: we could, for example, write the corre- 
lation matrix as (5(p a Sip b ) oc M ab , where M is a real 
positive symmetric matrix. One can then diagonalize M 
as M = UDU*. Redefining scalar field fluctuations as 
5(p ^ Sip = UStp will recover Eq. (|4]). 
Now, we impose the third condition: 

3. Truncate the SN expansion [Eq. j3j] at the or- 
der of Scp\ i.e., C = NJipt + \N ab 8<p a Jipl + 
j,N abc S^ b J^ + i,N abcd 5 V 'Z5v b Jv c Jvi Thus, 
we shall ignore the contributions in the power 
spectrum, bispectrum, or trispectrum coming from 
0(5ip 5 ). 

The 4th-order term is needed when we calculate all of the 
1-loop contributions in /nl and tnl- In the following, we 



Here, the dots in Eq. (JSJ) include the higher-order terms 
such as N 2 bcc V* 2 . This is a nice way of writing the power 
spectrum etc., as the results do not include the ultra- 
violet cutoff, fcmax, explicitly: the cutoff can be absorbed 
by redefining the derivatives of N. 

As we can take L such that kL 3> 1, 'Pioop is essentially 
a constant factor, rescaling the overall amplitude of the 
power spectrum without destroying the observed scale 
invariance of the power spectrum. Without loss of gener- 
ality, we shall take k to be the usual normalization scale 
used by the WMAP collaboration, fc = 0.002 Mpc -1 . 

Kawakami et al. have derived the exp ressions for 
/nl and tnl up to the 4th order (also see [14|)- These ex- 
pressions are again simplified by using the redefinition of 
the derivatives of N and ignoring the higher-order terms: 
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N a N a + Tr(N 2 )V* Hk L) 

N a N b N ab + (Tr(N 3 ) + 2N a N bc N abc ^j V, ln(fc i) 

(10) 
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N a N ab N bc N c + (2N a N ab N cd N bcd + Tr(A> 4 ) 
2N a N bc N bd N acd + N a N b N acd N bcd )v* Hk L) 



M ab = N ac N cb + N c N cab , we find 
(Tr(7V 3 ) + N a N bc N abc ) 2 P 2 ln 2 (k L) 



(11) 



where N abc = Nab c +\N a h cd dV* ln(A; max L). Although the 
loop terms of the bispectrum and trispectrum have terms 
like ln(fcbi), \n(k t L) and ln(fc p L) where k b = minjfci} 
with i = {1,2,3} or {1,2,3,4} , k t = min{fc l , \kj + fc;|} 
with = {1,2,3,4} and ln(fc p L) ~ ln(fcji) ~ 

ln(|fcj- + ki\L) with = {1,2,3,4}, we assume that 

these are similar to ln(fco^), i-e., \n(koL) ~ \n(k b L) ~ 
ln(fc f i) ~ ln(fcpL). From now on, we shall remove the 
tildes from the equations, i.e., N — > N. 

Now, we are ready to derive the new inequality. First of 
all, we use the inequality between arbitrary real numbers 
a and (3: a 2 + (3 2 > i (a + (3) 2 . Choosing a and /3 as 



< 



(iV a A^ Q ) 4 (l + P loop ) 4 
(Tr(iV 4 ) + 2N ac N cb N dab N d + N c N cab N abd N d ) V, ln(k L) 



(N a N a r(l+V l0 ^ 



loop } 



(15) 



where we have also used Pioop/(l + Aoop) < !• Here, 
let us reconsider the effect of our approximation that all 
the logarithmic factors are similar: ln(fcoL) ~ h\(k b L) ~ 
ln(fc f i) ~ ln(fcpL). If we relax this assumption, then 
we should replace ln(/co-L) in the right hand side of 
Eq. ([15]) with ln(feo-C') -> \n(k t L)R, where R = 
\n 2 {k b L)/hi{ktL) ln(fcpL). Therefore, our approxima- 
tion is valid also when the geometric mean of ln(fc t L) 
and ln(kpL) is similar to \n(k b L) (but not necessarily 
\n(k t L) ~ ln(fc p L)). 

Collecting these results, we obtain 



a=[N a N a (l 
(3 = [N a N a (l 

we find 



Vl oop )] [N a N b N ab + N a N bc N abc V* \n(k L)) 

Vi oop )]- 2 [Tt(N 3 ) + N a N bc N abc ] V, ln(fcoi), 

(12) 



1 f 6 

r NL + (2 loop) > - ( — ,/nl 



(16) 



[AUV Q (l+P loop )r 



N a N b N ab + N a N bc N abc V* ln{k L) 



(13) 



+ (Tr(iV ) + N a N bc N abc ) V 2 \n 2 (k L) > - f -fi 



1 f§ 



/NL 



Next, pick up the first term of the LHS in (fTSj) . and 
use the Cauchy-Schwarz inequality. When we define the 
inner product of arbitrary vectors v a and u b as (v,u) = 
J2 a VaUa ' then the Cauchy-Schwarz inequality leads to 
(v,u) 2 < (v,v)(u,u). Choosing v a and u a as v a = N a 
and u a = N b N ba + N bc N abc P* ln(k L), we find 



where the "2 loop" term is the last term in the RHS of 
Eq. (|T4"]) . This result shows that, when we allow our- 
selves for completely general models in which this par- 
ticular 2-loop term can become important, the Suyama- 
Yamaguchi inequality, tnl > (f/nl) > m ay be violated 
badly. This illustrates the limitation of this inequality. 

Still, from a model-building point of view, it is reason- 
able to assume that the 2-loop terms are sub-dominant 
compared to the tree or 1-loop terms; otherwise, we 
would have to require fine-tunings between the deriva- 
tives of N. Let us then study the consequence of ig- 
noring this particular 2-loop term. We shall impose the 
following conditions: 



(N a N b N ab + N a N bc N abc V* ln(fc L) 
(AUV a ) 4 (l+P loop ) 4 



N ab N abc N cde N de r 2 ln 2 (fc L) 
N b N ba N ac N c 

N ab N abc N cde N de V 2 \n 2 {k L) 
N b N ba N ac N c 



«1, 



N d N da N abc N bc V* Hk L) 



N b N ba N ac N c 



(17) 



< 



NbNbaNgdNd + 2N d N da N abc N bc V, ln(k L) 
(N a N a ) 3 (1 + Vioopf 
N ab N abc N cde N de V 2 \n 2 (k Q L) 
(N a N a )Hl+V loop ) 3 ' 



(14) 



where we have also used 1/(1 + 'Pi 00 p) < 1 with "Pi oop > 
on the RHS. Note that the last term on the RHS is a 
2-loop contribution, which becomes important later. 

Finally, pick up the second term of the LHS in (fl"3"| . 
and use the Cauchy-Schwarz inequality again: for arbi- 
trary real symmetric matrices M, L, we have Tr 2 (LM) < 
Tr(M 2 )Tr(L 2 ). Choosing L and M as L ab = N ab and 



The first condition is (tree) ^> (2-loop), and the second is 
( 1-loop) > (2-loop) for the terms in the RHS of Eq. (fl4| . 
Interestingly, from the Cauchy-Schwarz inequality for 
N ab N b and N abc N bc , we find 

( N d N da N abc N bc P, In(Aoi) \ 2 



N b N ba N ac N c J 
N ab N abc N cde N de Vl ln 2 (fc i) 



< 



N b N ba N ac N c 
N d N da N abc N bc V, ln(k L) 



N b N ba N ac N c 



(18) 
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from which we obtain the following bound on a particular 
form of 1-loop contributions: 



N d N da N abc N bc V* Hk L) 



N b N ba N ac N c 



« 1. 



(19) 



As a result, if we ignore the last term in the RHS of 
Eq. (|14l) . we must also ignore the second term, leaving 
only the tree-level term in the RHS of Eq. (fT4"|) . This is a 
peculiar feature of these terms, whose physical meaning 
is not clear. 

In any case, provided that the following additional con- 
dition is met: 

4. The 2-loop contributions are sub-dominant com- 
pared to the tree-level or 1-loop contributions (or 
at least the particular 2-loop term in the RHS of 
Eq. (fT4"]) is small compared to the others), 

we finally arrive at the new inequality: 

2 



1 /6 



2 V 5 



r/i 



NL 



(20) 



which is the main result of this paper, and is valid as 
long as the 2-loop contributions are small. This result 
generalizes the Suyama-Yamaguchi inequality (which in- 



cluded only the tree- level terms) as well as Ref. [14 1 
(which included up to the second-order terms). This 
relation can have a logarithmic scale dependence via 
R = ln 2 (k b L)/[ln{k t L) hx(k p L)]. 

What are the implications for inflation? In principle, 
if the following relation is observed, 



1 

- < 



tnl 



{tJnlY 



< 1, 



(21) 



then it implies that there was a non-negligible contribu- 
tion from 1-loop terms, which would help constrain the 
physics of multi- field models (via the forms of N). How- 
ever, the most interesting case would be the observation 
of a complete violation of the inequality, i.e., 



1 /6 , 
tnl "C - I ^jnl 



(22) 



which implies that inflation cannot be responsible for 
generating the observed fluctuations, provided that (1) 
scalar fields are the source of fluctuations; (2) fluctu- 
ations at the horizon crossing are scale invariant and 
Gaussian; (3) the evolution of fluctuations obeys the 6N 
formalism; and (4) the 2-loop contributions are small. 

We may not be so far away from testing this predic- 
tion. If the value of /nl is as large as what is implied 
from the current data, /nl ~ 30, then the threshold 
value, tnl ~ 650, is close to the 2- a limit expected from 
Planck 2l|, |22j. The large-scale structure observations 



should also help improving the limits on tnl [HI . There- 
fore, in the event that Planck sees /nl (thus ruling out 
single-field models [24]), one of the two things can hap- 
pen: (1) t N l is also detected in excess of tj(§/nl) 2 , con- 
firming predictions from multi-field models, or (2) tnl is 
either not detected, or detected below ^(|/nl) 2 , ruling 
out most of the multi-field models that satisfy the above 
4 conditions. This argument [llj and our result provide a 
strong science case for measuring the local-form trispec- 
trum of the cosmic microwave background as well as that 
of the large-scale structure of the universe. 
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